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Abstract. Dust Dynamics Simulations utilizing a dynamical screening approach are 
performed to study the effect of ion-streaming on the self-organized structures in a 
three-dimensional spherically confined complex (dusty) plasma. Varying the Mach 
number M - the ratio of ion drift velocity to the sound velocity, the simulations 
reproduce the experimentally observed cluster configurations in the two limiting cases: 
at M = strongly correlated crystalline structures consisting of nested spherical shells 
(Yukawa balls) and, for M > 1, flow-aligned dust chains, respectively. In addition, 
our simulations reveal a discontinuous transition between these two limits. It is found 
that already a moderate ion drift velocity (M w 0.1) destabilizes the highly ordered 
Yukawa balls and initiates an abrupt melting transition. The critical value of M is 
found to be independent of the cluster size. 



Ion- Streaming Induced Order Transition 



2 



1. Introduction 

Collective dynamical effects such as wakefields are omnipresent in nature. One of the 
most common phenomena is the V-shaped 'Kelvin wake pattern' created by watercrafts 
moving through water. In a plasma, similar wake pattern give rise to an effective 
attractive interaction between particles that otherwise mutually repel each other. In 
recent years, plasma wake fields which are driven by an ultrashort laser pulse are used 
for the acceleration of electron beams where it is expected that rates that are more 
than a thousand times higher than those achieved in conventional accelerators may 
be achievable [U [2] • In an entirely different context - conventional superconductors - 
similar attractive forces are created by the dynamically polarised background medium 
giving rise to the Cooper pairing of electrons. 

An excellent test system for many-particle correlations and structure formation 
with pronounced streaming and wake effects are complex (dusty) plasmas [3J. In the 
plasma bulk of an radio-frequency (RF) discharge, where the plasma flow is negligible, 
spherical 3D plasma crystals (Yukawa balls) with a nested shell structure have been 
created [1]. A completely different type of dust structure is observed in the regime 
of supersonic ion flow, M > 1. In the plasma sheath region, where strong electric 
fields are present, the grains form vertically aligned 'dust molecules' and, with the 
help of a confining potential, long particle strings can be created, e.g. [5], El [7]. The 
vertical alignment of particles is not explainable with a purely repulsive grain interaction. 
Therefore, it was recognized early that the highly charged grains have a focusing effect 
on the streaming ions and a positive ion space charge is accumulated in the wake 
below the grains [S]. Experimentally, the existence of an attractive wake potential 
and a resulting non-reciprocal grain interaction was demonstrated by means of laser 
manipulation experiments on particle strings [5] and dust molecules [6l[9j. 

While the two limiting cases (Yukawa balls and flow-aligned dust strings) have 
been explored experimentally and theoretically in great detail, see e.g. |3], [10] and 
references therein, a key question has- to our knowledge -not been addressed so far: 
how does the structural order transition proceed between both limits when the ion 
flow is gradually increased and wake effects increasingly take charge? Experimentally, 
this question may be hard to answer since the ion-streaming velocity is neither directly 
accessible nor separately tunable without changing other plasma parameters at the same 
time. Therefore, computer simulations may help to answer these questions. However, 
a self-consistent numerical simulation of a nonideal multi-component plasma in non- 
equilibrium is highly challenging due to the drastic differences in the time scales of the 
plasma constituents: electrons, positively charged ions, neutral atoms, and negatively 
charged 'dust' grains of micrometer size. Utilizing the 'Dynamical Screening Approach' 
(DSA) this multiscale problem can, however, be effectively reduced by an approximate 
(statistical) description of the lighter plasma constituents [H HH [121 [131 [IH [Ej • The 
DSA allows for an accurate description of essential equilibrium and non-equilibrium 
plasma properties including plasma streaming induced wakefield oscillations, ion and 
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electron thermal effects as well as coUisional and Landau damping. On this basis, the 
dust dynamics will be studied by first principle Langevin Dynamics simulations, which 
allow for a full dynamic treatment of the dust pair interaction. 

This article is outlined as follows: In section [2] the DSA is employed to incorporate 
the dynamical screening properties of the streaming plasma background into an effective 
electrostatic potential of a single grain. In section [3] we then use these results for many- 
particle dust dynamical simulations which reveal the collective particle behavior for a 
wide range of ion-streaming velocities. A final discussion of the ion-streaming induced 
order transition follows in section HI 

2. Dynamically Screened Grain Potential 

Collective dynamical effects of the highly charged, massive dust particles are mainly 
determined by their electrostatic interactions. For the computation of the electrostatic 
potential of a single grain we consider plasma parameters that match the typical 
experimental conditions: plasma density rig = rii = 10^/ cm'^, electron temperature 
Tg = 2.5 eV, ion temperature Tj = 0.03 eV. A moderate Argon gas pressure of 
p = 15Pa corresponds to an ion-neutral collision frequency of z/j„ = 0.2 Up. [16], where 
Wp. = \J Tiiof / {e^mi) = 2.1 ■ 10^ Hz is the plasma frequency of singly charged Argon ions 
of mass mi = 6.634 ■ 10~^^kg. 

Assuming a weak (linear) response of the plasma to the presence of the dust grain, 
the dynamically screened grain potential in a streaming plasma is obtained in real space 
by means of a 3D Fourier transformation 

where r denotes the distance from the grain's center of mass [T7j. For the calculation 
of the wake potential we use a grain charge of = — IO^Cq, but the potential can 
be rescaled linearly to any other charge value, see section [3j The ambient plasma is 
represented by a longitudinal dielectric function which includes a directed ion flow Uj 
and Bhatnagar-Gross-Krook (BGK)-type ion-neutral collisions [H] 

Here, we employ the substitution Q = (k • Uj -f- iuin) / {\/2kvTi) the plasma 
dispersion function Z(z) = it~^ dt '^^^i"* fT^. Further, we introduced the Debye 
length of electrons Ad,. = 1175.41 /im and ions Ad, = 128.76 /xm, respectively, which 
are defined as Ad,^ = (sokBTa/riaqaY^'^ ■ Note that screening by the electrons is 
purely static, since their thermal speed exceeds their field-induced drift. The ion 
thermal speed = v^/^B^i/mj = 270 m/s is well below the ion sound (Bohm) speed 
Cs = ^Jk-QTe/rrii = 2460 m/s. The ion streaming velocity is expressed in terms of the 
Mach number M = |uj|/cs. 



1 + 0^(0) 

1 + ^^(c^: 



(2) 
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Figure 1. (Color online) Dynamically screened dust (wake) potential "I>(r) = z) 
as obtained from Eq. ([T]), = 0.128 mm. The grain is located at {p, 2:} = {0,0}. 
The ions stream along the z-direction with M = 0.01,0.08,0.16,0.33,0.66, and 1.0, 
respectively, giving rise to a wake structure behind the grain. The values of the 
potential along the z-axis are shown also on Fig. [2] 

In Figs. [1] and [2| we present the dynamically screened Coulomb potential for 
a single grain in the relevant range of ion flow velocities Uj. When ion drift and 
collisions are negligible, that is M — and z/j„ — )■ 0, Eq. ^ reduces to static 
electron and ion shielding of the Coulomb potential and the grain potential becomes 
the isotropic Debye-Hiickel (or Yukawa) potential $Yuk(^) = "Jd e~''/^°/(47reo^), where 
Ad = (Aj3^ + Aj3^^)^^/^ = 128.00 /xm is the combined Debye length of electrons and ions. 
Fig. [1^) shows a 3D plot of the electrostatic potential $(r) = $(r, z) for M and 
finite pressure. 

A stationary flow of ions along the z-direction leads, however, to strong deviations 
from the Yukawa potential by giving rise to an anisotropic, oscillating wake structure 
behind each grain, as shown in Fig. [T|D)-f). With increasing M the amplitude of the 
ion focus increases and reaches its maximum for M f» 0.5, where $max = 43 mV. This 
means that essential wake effects are already present well below the supersonic regime. 
A consideration of an increasing dust charge qd with M does not change the wakefield 
significantly, see OML results in Ref. [T3]. In Fig.jlji), the first minimum and the second 
maximum of the potential become clearly visible in the wake structure. As a general 
trend, it is found that with increasing M the primary ion focus and all higher extrema 
shift linearly away from the grain. In situations, where the streaming velocity is close 
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Figure 2. (Color online) Electrostatic grain potential <i>(r, z) in streaming direction z 
at the radial position p = x = y = [cut through Fig. [ij. 



to or exceeds the sound speed, the ion focus flattens out and a more and more stretched 
Mach cone forms. More details on the numerical implementation, the wake structure, 
and a critical assessment of the DSA will be published separately in Ref. [U] 

3. Dust Dynamics Simulations 

The dust grains are embedded in a partially ionized plasma and exposed to collisions 
with the neutral gas atoms. Instead of treating the dust-neutral collisions explicitly, the 
Langevin scheme provides a stochastic approach, which involves the neutral g 
large heat bath for the dust particles. The equation of motion for the k-th dust grain 
reads 

mdTk = -W^^{r, t) - uliTLdYk - ydnmdVk + ffc(t) • (3) 

Since Maxwell's equations are linear, the total dynamic potential acting on the k-th 
grain is just the superposition of the electrostatic (wake) potentials created by all other 
{Nd — 1) grains 

Vf{v,t) = Y,<lMr.t) ■ (4) 

A screening independent parabolic confinement is explicitly included in Eq. (|3]). It 
mimics the combined effect of gravity, thermophoresis, and electric fields which in the 
experiments with Yukawa balls are used to achieve an approximately isotropic harmonic 
confinement for the dust grains. [TH] The friction term emulates a damping mechanism 
and a stochastic force, in(t), accounts for random heating effects. The random force and 
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Figure 3. (Color online) Sim- 
ulation snapshots of the (sta- 
tionary) state of Nd = 46 par- 
ticles and the ion-streaming 
velocities M = 0.01,0.08,0.16, 
and 1.00. The plasma streams 
from top to bottom. The color 
code for M = 0.01 and 0.08 
denotes the particle's radial 
distance from the trap cen- 
ter. Approximate cluster size 
(in z-direction) : a) 2.2 mm, 
b) 3.0 mm, c) 6.3 mm, and d) 
7.0 mm. 
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the (neutral gas) friction term are related by the fluctuation-dissipation theorem 

(4"(t)ff (t')) = 2m,u,,,kBTJkiSap6{t - t') , (5) 

where a, (3 G {x,y,z} and k,l G {l,...,Nd}. The standard Langevin algorithm 
employed here is extensively described in Ref. [19]. Connecting Langevin dynamics 
with the dynamical screening approach, the DDS scheme allows for inclusion of flnite 
temperature effects as well as a full dynamical treatment of all dust-dust interactions on 
the basis of the effective wake potential Eq. ([T]) [THl . Besides the plasma parameters 
used for the computation of the grain potential, additionally the mass = 9.1-10"^^ kg, 
charge = — 6000eo, and radius Rd = 2.43 yum of A^^^^ identical grains as well as the 
trap frequency uq = 7.0 Hz enter the simulation as input parameters. A dust-neutral 
collision frequency of i/dn = 5.3 • \/mn/kBTn = 19.1 Hz drives the dust temperature 
towards the neutral gas temperature Tn = Ti. 

In the following we will investigate the effect of the only remaining free parameter 
in the considered model system: the relative streaming velocity of the ions with respect 
to the dust in terms of the Mach number M. In the limiting case of no ion flow, i.e. 
M = 0, the grain potential is of Yukawa-type. Starting from a random initial state, 
the grains equilibrate during the simulation due to frictional drag against the stationary 
neutral gas. Finally, the sub-system of dust grains reaches a strongly coupled state, 
where the dust forms a spherical Yukawa ball with a nested shell structure. Fig. [3^) 
shows the cluster with Nd = 46 grains which exhibits a meta stable (1, 12, 33) shell 
conflguration for M = 0. The highly symmetric icosahedral 12 particle conflguration in 
the inner shell is known to be particularly stable and was found to be energetically most 
favorable for isotropically screened Coulomb particles in a parabolic conflnement. |21] . 
Earlier simulations in the frame of the one-component plasma model with a statically 
screened Coulomb potential ^vukl'") have proven to yield very good agreement with the 




Figure 4. (Color online) Pair 
distributions for iV = 18 grains 
at the streaming velocities M = 
0.01,0.08,0.16, and 1.00. 



Figure 5. (Color online) 
Vertical (string) order parameter 
as defined in Eq.Q for different 
particle numbers N. 



experiments in view of the structure [221 [23], dynamics [3], and melting behavior [21] of 
finite dust crystals. 

Next we consider the case of supersonic ion-streaming. When the ions are streaming 
along z-direction from top to bottom, they are deflected by the highly charged dust 
grains and induce an ion focus directly below the grain. Other grains downstream 
are attracted to this positive space charge. As shown in Fig. [sji), the ion wakefield 
attraction leads in the case of = 46 to the formation of four vertical particle 
chains. Because of the asymmetric nature of the (non-conservative) wake potential 
along the z-axis, the interaction between an upstream and a downstream grain is 
non-reciprocal [H] [I3]. Kinetic energy of the streaming ions is coupled into the non- 
Hamiltonian dust system, which leads to self-excited vibrations and eventually to the 
instability of long chains [131 [23 ES]- This effect is indicated in Fig. |3]i), where the dust 
particles in the upper part of the particle strings are highly ordered, but in the lower part 
grain dislocations due to a downward increasing, streaming-induced dust temperature 
are observed. 

While in the two limits of M = and M > 1 the grains are strongly coupled, 
the transition between both limits is found not to be necessarily smooth. Instead, a 
relatively weak ion flow of the order of M = 0.1 is found to completely destabilize 
the highly ordered dusty plasma state and the grains become spatially delocalized, 
see Fig|3]D). At M = 0.16, the formation of a short-range ordered 'doublet state' is 
observed. The vertically aligned pairs of grains move through a weakly coupled, liquid- 
like state. As observed e.g. in laser manipulation experiments [HI [6], the downstream 
particle strictly follows the upper particle's motion. In contrast, if the lower dust particle 
experiences a lateral kick by other grains, the (vertical) center axis of the dust molecule 
is tilted sideways. If the lower grain reaches a sufficient height with respect to the 
upper grain such that it leaves the attractive wakefield, the laterally repulsive grain 
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interaction pushes the grains in the opposite directions and the binding between the 
grains is destroyed. An ion flow of M = 0.33 is just above the threshold, where 'triplet 
states' become stable. The lowest grain shows a strongly oscillatory behavior and other 
free grains or complexes are too inert to follow that motion. Moreover, the (electrostatic) 
interference of other grains easily destroys the binding of the lowest grain. With further 
increase of M more and more grains line up behind each other and form longer, streaming 
aligned particle chains. 

A similar ion-streaming mediated order transition is found for clusters of different 
size. The time-averaged pair distribution function for A^^^^ = 18 grains in Fig. |4] gives a 
qualitative picture of the ion-focus induced melting transition and subsequent structure 
built-up upon changing M. The distinct peak structure in the two limits M = 0.01 and 
M = 1.0 is a clear signature of strong correlations in the underlying dust cluster. That 
is a Yukawa ball with a (1, 17) configuration in the static case. At ion Bohm speed, 
the confining potential makes it possible to have (two) nine-particles long fiow-aligned 
strings. These strings exhibit self-excited oscillations which grow to a large amplitude 
at the bottom end of the strings, as described in Ref . jT3] . These ion-streaming induced 
oscillations are a reason for the broadening of the pair distribution at large grain 
distances. The pair distribution function for M = 0.08 is (almost) structureless which 
indicates a fiuid-like state. For M = 0.16 we find a single sharp peak which signals the 
formation of temporal stable dust doublets. 

To quantify the character of the string formation, we make use of the Vertical Order 
Parameter ( VOP) which was recently introduced by Killer et al. [7] . The VOP captures 
the fraction of vertically aligned bonds with (p < 10° to the total number of bonds 



where (j) denotes the angle between the bond axis of two grains and the vertical z-axis. 
The VOP is equal to one, if the particles create a single chain of particles (which may 
exhibit transverse oscillations [25]). It is zero, when there is no vertical order such as in 
a fiat 2D system. The (time-averaged) VOP in dependence on M is plotted in Fig. [sjfor 
different particle numbers N^. As a general trend, none of the spherical clusters observed 
in the static case exhibits a vertical particle alignment. This finding is in agreement 
with Ref. [7]. Independent of the exact particle number and finite size effects, the 
VOP shows a sharp jump around c<j/10, when the ion flow increases from M = 0.08 to 
M = 0.16. This crucial effect is due to the formation of spatially correlated doublet 
states. When M is further increased, the formation of triplets and subsequently the 
formation of long particle chains is observed. Vertical alignment is most pronounced for 
M = 1.0. For M = 1.33 the value of the VOP is slightly reduced due to a conflgurational 
rearrangement: going from M = 1.00 to M = 1.33, there are for = 46 [N^ = 18) 
seven (three) instead of four (two) vertical strings. 

At this point, the question concerning the mechanism for the ion-streaming induced 
melting process remains. In the static case, the mutual (isotropic) grain repulsion is 



5, 



tot — 



N,{N, - l)/2 



VOP 




(6) 
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compensated by the external confining potential. That stabilizing interplay allows for 
strong correlations. At M = 0.08, however, the grain potential in downstream direction 
becomes non-monotonous and relatively flat. Consequently, the repelling force between 
the grains vanishes and the mechanism for the correlation built-up breaks down. Only 
with further increase of the ion flow, a different kind of order process emerges which 
does not stabilize the system as a whole, but rather individual grains agglomerate to 
larger complexes (doublets, triplets, etc.) due to wakefield attraction. That is, a global 
stabilization process is replaced by an local one. 

4. Summary and Conclusions 

The influence of an external electric field on the dynamical alignment (laning) of like- 
charged particles is a topic of high current interest, e.g. in the physics of colloids [2Z]. 
In dusty plasmas physics, however, the important effect of the (field-induced) ion-flow 
on the non-equilibrium structure formation has not been systematically analysed for 
the full relevant range of Mach numbers M so far. To this end, we have performed 
first principle Langevin-type Dust Dynamics Simulations with a dynamically-screened 
Coulomb potential obtained from linear response theory. 

The model considered here is found to reproduce the experimentally observed 
particle arrangements, that is Yukawa balls in the static case (M = 0) and flow-aligned 
particle chains in the supersonic regime (M > 1). To focus on the explicit effect of a 
finite Mach number M, the influence of M on other plasma parameters such as the grain 
charge or trap frequency was neglected. The simulations reveal that even a relatively 
moderate ion-flow initiates an abrupt melting of the highly ordered (equilibrium) many- 
particle state. The observed, abrupt melting process resembles a non-equilibrium phase 
transition. With further increase of the ion flow, the delocalized grains start to form 
flow-aligned complexes: doublets, triplets, and longer particle strings. For M > 1, a 
re-entrance into a strongly coupled, non-equilibrium many-particle state consisting of 
several (oscillating) chains is observed. A quantitative analysis of this non-equilibrium 
phenomenon by means of the Vertical Order Parameter shows that the order process 
appears in the same M-range for clusters of different size. Since the underlying, wake- 
induced, order mechanism is of local nature, we expect a similar ion-focus induced order 
transition to appear in similar form also in unconfined macroscopic systems. 

Since dusty plasmas are often considered as a test system for many-particle 
correlations, the finding that even weak streaming can trigger unexpected many-particle 
behavior, may be of relevance also for other fields of physics. For instance, in warm 
dense matter physics, electron streaming effects are often neglected in the theoretical 
description. Using a quantum (e.g, Mermin) dielectric function, the DDS scheme can 
be directly applied also to quantum systems [20] . 
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